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Abstract 

We report on an exact method for global fits of the CKM matrix by using 
the necessary and sufficient conditions the data have to satisfy in order to 
find a unitary matrix compatible with them, and this method can be ap- 
plied to both quark and lepton CKM matrices. The key condition writes 
— f < cos 5 < 1 where S is the phase that encodes the CP violation, and it 
is obtained when one describes the CKM matrix in terms of four moduli 
of its entries. We use it to reconstruct many approximate unitary matri- 
ces from the PDG data, and by using the double stochasticity property 
of the unitary matrices we find a strong dependence of sin 2/3 on |f7 u b|, 
that implies also strong correlations between \Uij | themselves, proving the 
constraining power of this approach. 



1 Introduction 

Recently we have introduced a unitarity condition method for constraining 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix entries that parameterize the 
weak charged current interactions of quarks in the Standard Model. The de- 
termination of the CKM matrix is a lively subject in particle physics and the 
determination of its four independent parameters that govern all flavor changing 
transitions of quarks and leptons in the Standard Model is an important task 
for both experimenters and theorists. The large interest in the subject is also 
reflected by the workshops organized in the last years whose main subject was 
the CKM matrix 00 HE]. 

Although it seems that at the level of experts there exists a consensus con- 
cerning the method of extracting from experimental data information about the 
CP violating phase and the invariant angles of the unitarity triangles 0-^11], the 
usual method of considering only the orthogonality of the first and third columns 
has some weak points as it was shown in pQ. Here we give more arguments in 
favor of our novel approach; the main point of the unitarity method consists 
in taking as independent parameters four moduli of the CKM matrix entries, 
and finding from them constraints on the phase S. An alternative approach 
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could be the use of four independent invariant phases or any combination 
of moduli and phases, the sole condition to be fulfilled being their indepen- 
dence. In fact the standard unitarity triangle method 0-^01 as it is used in the 
present days allows the determination of only two independent invariant phases, 
showing that using only one orthogonality relation is not enough for a complete 
determination of the CKM matrix. Moreover we have constructed in P a toy 
model matrix that satisfies exactly what we call the double stochasticity prop- 
erty, see below, its first and third columns arc orthogonal, and the matrix is not 
unitary. Until now the toy matrix was considered as an academical curiosity, 
but in fact it stresses the necessity to consider at least another unitarity triangle 
for the determination of others two phases, and a good candidate could be that 
determined by the orthogonality of the first and the third rows. 

Our approach provides the necessary and sufficient condition the data have 
to satisfy in order to find a unitary matrix compatible with them; this condition 
is given by — 1 < cos<5 < 1, where 5 is the phase entering the CKM matrix. In 
section 2 we construct our theoretical model and give a few numerical examples 
that show the constraining power of our approach. In section 3 we make a 
comparison between the standard unitarity triangle approach and ours and do 
some Monte Carlo simulations stressing the numerical problems that have to be 
solved when doing global fits. In section 4 we test the unitarity property of the 
experimental data by using only the information provided by the Particle Data 
Group (PDG) 22]- mill and we found that there are many unitary matrices 
leading to a phase S close to tt/2. We use the double stochasticity property 
and the unitary matrices found in the fits to construct continuous families of 
unitary matrices, and as an unexpected result we uncover an explicit dependence 
of sin 2/3 on | U u b \ . Thus we can use the experimental data on sin 2/3, see e.g. [H] . 
and the mean unitary matrices obtained from Hagiwara et al. data |12j . and 
respectively Edclman ct al. data to determine unitary matrices compatible 
with all these data, by reversing in some sense the usual way of getting bounds 
for sin 2/3 from fits by using the moduli \Uij\. Working in this way all the 
found matrices as well as their ±lcr companions are unitary. The paper ends 
by Conclusions. 



2 Theoretical model 

We use the standard parameterization advocated by PDG that we write it in a 
little different form 

(C12C13 C13S12 si 3 \ 

-C 2 3Sl2e~[ S - C12S23S13 Ci2C2 3 e" 1 ' 5 - S12S23S13 S23C13 (1) 
Sl2S23e~ i<5 - C12C23S13 -C 12 S 2 3e~ lS - S12C23S13 C23C13 / 

with dj = cos 9ij and = sinfl^ for the generation labels ij = 12, 13, 23, and 
S is the phase that encodes the breaking of the CP-invariance. The rc-phasc 
invariance is equivalent with choosing the phases of five arbitrary CKM matrix 
entries equal to zero, and we used this property for rewriting the standard 
parameterization |15| in the above simplest form. This form depends on four 
phases, those of U21, U22, t/31 and U32, instead of five as in the usual approach. 
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From experiments one measures a matrix whose entries are positive 



V 2 d V 2 V 2 (2) 

v t \ v 2 V 2 J 

entries that, in principle, can be determined from the weak decays of the relevant 
quarks, and/or from deep inelastic neutrino scattering As the notation 

suggests we make a clear distinction between the unitary CKM matrix U and 
the positive entries matrix V provided by the data. Of course the data are in 
general given in terms of some functions fk{Vij), k = 1, . . . , N, that depend on 
the V entries. 

The main theoretical problem is to see if from a matrix as one can 
reconstruct a unitary matrix as For that we need a relationship between 
the theoretical object and the experimental data @, and it is 

V = \U\ 2 

that leads to the following relations 



* ud 


- r 2 r 2 V 2 

— c 12 c 13' v us 


_ s 2 „2 y2 _ 2 
— 6 12 c 13' v ub ~ 6 13 


Vcb 


s 2 r 2 V 2 


2 2 
— c 13 C 23i 


v 2 d 


2 2,2 
— s 12 c 23 "+" s 13 


s 23 c 12 + 2S12S13S23C12C23 COS S, 


V 2 

cs 


2 2,2 
— C 12 C 23 ' S 12 


s 13 s 23 ~ 2S12S13S23C12C23 COS S, 


1 Id 


2 2 2, 
^ S 13 C 12 C 23 "r 


S 12 S 23 _ 2S12S13S23C12C23 COS (5, 


v t 2 


2 2 2, 
— S 12 S 13 C 23 + 


c 12 s 23 + 2S12S13S23C12C23 COS S 



From the relations 10 one gets 

V^-1 = 0, j = u,c,t 

i—d,s,b 

E V^-1 = 0, j = d,s,b (4) 

i—u,c,t 

that is the weakest form of unitarity. We stress that the above relations does 
not test the unitarity, as it is usually stated in many papers; they are necessary 
but not sufficient conditions. The class of positive matrices satisfying Eqs.Q is 
considerable larger than the class of positive matrices coming from unitary ma- 
trices. The set (@J is known in the mathematical literature as doubly stochastic 
matrices, and the subset coming from unitary matrices V 2 - = \Uij\ 2 is known 
as unistochastic ones [EJ. The double stochastic matrices have an important 
property, they are a convex set, i.e. if V\ and V2 are doubly stochastic so is their 
convex combination a V\ + (1 — a)V2, a G [0, 1], as it is easily checked. 

The first problem to solve is to find a necessary and sufficient criterion for 
discrimination between the two sets. The above relations provides the necessary 
and sufficient conditions the data have to satisfy in order the matrix (0 comes 
from a unitary matrix, and the most constraining condition is 

-l<cos<S<± (5) 
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Since in relations © d enters only in the cosine function we can take d £ [0, tt] 
without loss of generality. Using the central numerical values from PDG data 
or from the published fits 0-0, one finds values for cos S outside the 
physical range [—1,1], or even complex as it was shown in pQ. The last four 
relations J3J provide us formulas for cos 8 and these formulas have to give the 
same number when comparing theory with experiment, by supposing the data 
come from a unitary matrix. Their explicit form depends on the independent 
four parameters we choose to parameterize the data, and we will always choose 
these parameters as four experimentally measurable quantities, i.e. V 2 y In 
fact there are 57 independent groups of four independent moduli that lead to 
165 different expressions for cosi5. Depending on the explicit choice of the four 
independent parameters we get one, two, three or four different expressions for 
cos 5; e.g. if we take V u d, V us , V cd , V cs , as independent parameters we get 



V ud 1 us V ud 1 us 

and from the sixth Eqs.10 we have 



V 2 , + V 2 ,V 2 , + V 2 V\ + V\, - V 2 V\ 

y ud 1 y cd y ud 1 cs ud 1 v ud y cs y ud 

+V 2 -V 2 jV 2 -V 2 V 2 +V 2 jV 2 jV 2 -V a ~~V 2 V 2 ^V 2 -{-VlV 4 (7) 

c _ 1 us y cd us y cs v us 1 y cd ud us y us v cs v ud ns 1 v cd us \ ' J 

2v ud v us Vi-VS-K 2 . Vi-v 2 d -vi ^v- d + vi + v- d + v^-i 

In this case, other two independent formulas arc given by the last two equa- 
tions in relations J2Jl, and they have to give (almost) identical numerical results 
if the data are compatible with unitarity. In general the data will give different 
numerical values for the three functions expressing cos<5. If the independent 
parameters are V u d,V u b,V c d,V c b, i.e. we use the information contained in the 
first and the third columns as does the standard unitarity triangle, we obtain 
four different expressions for cos 6. In this case the mixing angles 6ij are given 

by 



and we obtain for cos 5 the following formula 



cos S = 

V cd YMk - + V cb V ub + VclVl - VI + VMM + V ub V us + VclKlVl 

2 Vus V ub V cb VI - Vj b - VI VI - V 2 b - VI 

Looking at equations we see that the expressions defining the mixing 

angles and phase 5 are quite different. Thus if the data are compatible to the 
existence of a unitary matrix these angles Sij and phases 6^' have to be equal, 
and these are the most general necessary conditions for unitarity; they can be 
written as 

< 4 n) < 1, 4? ] = s ii ] > m,n = 1,..., 57, cos^ (?;) = cos5 (i) , i,j = l,..., 165 
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The above relations are also satisfied by the double stochastic matrices, and the 
condition that separates the unitary matrices from the double stochastic ones 
is given by the relation i.e. —1 < cos<r 1 ' < 1. 

As a warning, what we said before can be summarized as follows: the unitar- 
ity property is a property of all the CKM matrix elements and not the property 
of a row and/or a column, and it implies strong correlations between Uij. 

To better understand the above considerations, let us consider the most 
general exact form of a double stochastic matrix that is obtained when we 
assume that the relations (@J are exactly satisfied, or in other words the V entries 
are known with infinite precision. Because almost all the unitarity triangle fits 
use U us , U u b, and U c b as input parameters we make the following notation 



\U us \=a, \U ub \=b, and \U, 



cb 



Since we need four moduli we choose as the fourth one, \U c d\ = d, and with they 
we form the matrix 



Vl — a 2 — b 2 a b 

s = \ d yi - c 2 - d 2 __c ] (io) 

Va 2 + b 2 -d 2 V-a 2 + c 2 +d 2 yjl - b 2 - c 2 

The squared matrix S 2 is an exact doubly stochastic matrix as it is easily 
checked. If we make the substitution V%j — > SV,-, all the 165 cos 5 expressions 
have the same form, namely 

r d 2 (l-& 2 ) 2 -a 2 (l-6 2 -c 2 )-& 2 c 2 (l-a 2 -& 2 ) , , 

COSO = ; ; (11) 

2abcVl -a 2 - b 2 Vl -b 2 -c 2 v ; 

From the condition 

cos 2 5 < 1 



we obtain the equivalent relation 

a 4 - 2a 4 c 2 - 2aW + a 4 c 4 + 2a 2 6 2 c 4 + 6 4 c 4 - 2a 2 d 2 + 2a 2 b 2 d 2 + 2a 2 c 2 d 2 
-2b 2 c 2 d 2 + 2a 2 b 2 c 2 d 2 + 2b i c 2 d 2 + d 4 - 2b 2 d± + b i d i < 



(12) 



that describes the physically admissible region in the 4-dimcnsional space gen- 
erated by the moduli a,b,c,d. For a similar result see |17| . where the authors 
proposed, for the first time, a parametrization of the CKM matrix in terms of 
four independent squared moduli, and they looked for constraints implied by 
unitarity. This result proves that four independent moduli do not determine 
a unitary matrix, even when the moduli satisfy the relations 10} ; this happens 
then and only then, when the relation JS} is satisfied, and the above relation 
separates the unitary matrices from the doubly stochastic ones in our choice of 
independent moduli. If we select other four moduli we find another form of the 
region. To see that the above constraint is very strong, let us take the median 
values from \1?>\ for a, 6, c, and then we find that cos 5 takes physical values if 
and only if 

0.223659 <d< 0.223958 

i.e. a very small interval from the allowed range 0.221 < d < 0.227. If we use 
instead of the central values for the above parameters the recommended values 
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given in |13| and write them as rational numbers, i.e. a = 0.22 = 11/50, b = 
0.00367 = 367/100000, c = 0.0413 = 413/10000, d = 0.224 = 28/125, one finds 



cos S = 



1160896633137657437088689 



28.6448 



4158110000000^94995845323482088321 



If we take a = 0.224, i.e. we modify it only with the small quantity 4 x 10 -3 
we find cos S ss 1.28, that shows that the fulfillment of unitarity imposes very 
strong correlations between all the entries of the CKM matrix. 

Now we define a test function that has to take into account the double 
stochasticity property expressed by the conditions Q and the fact that in gen- 
eral the numerical values of data are such that cos S depends on the choice of the 
four independent parameters and could take values outside the physical range. 
Our proposal is 



This proposal is our main contribution to the existing methods for fitting the 
CKM matrix entries, and it expresses the full content of unitarity. We stress 
again that the both conditions have to be fulfilled, x\ nas to take small values 
and cos 5 values have to be physically acceptable, since the relation Xi = holds 
true for both double stochastic and unitary matrices. The second component 
that takes into account the experimental data has the form 



where V%j is the numerical matrix that describes the experimental data, and a 
is the matrix of errors associated to Vy . In the last sum we use only the data 
coming from the first two rows because the entries of the third row arc not yet 
measured. The above form can be easily modified if we have other experimental 
information for some functions depending on Vij . 

The above expressions can be used to test globally the unitarity property of 
the experimental data and of the published fits. 

3 Comparison of the two methods 

Unitarity triangle approach exploits other consequences of the unitarity property 
of the matrix Q, namely the orthogonality relations of rows, and, respectively, 
columns of a unitary matrix. Although there are six such relations usually one 
considers only the orthogonality of the first and the third columns of U, relation 
that is written as 




(13) 




(14) 



U ud U* ub + U cd U* cb + U td U* tb = 



(15) 
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The above relation is scaled by dividing it through the middle term such that 
the length of one side is 1. The other sides have the lengths 



Ru 

Ri 



U udU* b 

Uc d U* cb 

u td u t * b 



b^\-a 2 -b 2 
cd 



U cd U* 



Va 2 + b 2 -d 2 Vl-b 2 - c 2 



(16) 



where we have written on the right hand side the R- values in our choice of the 
four independent parameters that lead to the double stochastic matrix i|lU|) • If as 
in the preceding case we compute the expressions on the right hand side using the 
values recommended in p|] we find R u = 0.387923, and R t = 4.53416i. If R u is 
in the "normal" range, Rt gets imaginary. Please remark that both approaches 
send the same signal, although in different ways: something is wrong with the 
recommended values. Using the above expressions one sees that Rt gets positive 
only if a 2 + b 2 - d 2 > 0. 

To have a quantitative estimate of the constraining power of both approaches 
we did a Monte Carlo simulation of the expressions appearing in and l|16|) 
using as input a = 0.2200 ± 0.0026, ,b = (3.67 ± 0.47) x 10~ 3 , c = (41.3 ± 
1.5) x 10~ 3 , d = 0.224 ± 0.012 and the results are shown in Fig[l]and FigEl 
The results show that only a tiny fraction of simulated data arc in the physical 
region. 




Figure 1: Number of Events as function of cos<5. The input values are: a = 
0.2200 ± 0.0026, b = (3.67 ± 0.47) x 10~ 3 , c = (41.3 ± 1.5) x 10" 3 , d = 0.224 ± 
0.012. On the left panel one sees the number of events in the physical region, 
cos S £ (1, 1). On the right panel are all the events and the interval of variation 
is cos S £ (—350, 450). The number of physical events is f=a 1% of the total 
simulated data, and there is a peak around cos 5 ~ 0. 
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2250 - 




R„ R, 2 



Figure 2: Number of Events as function of R u , and i?| . The input values are: 
a = 0.2200 ± 0.0026, b = (3.67 ± 0.47) x 10" 3 , c = (41.3 ± 1.5) x 10" 3 , d = 
0.224 ± 0.012. R u takes values within the interval R u G (0.2, 0.6), and R 2 t G 
(—250, 250), and only 37.5% of them are positive. If we impose the triangle 
inequality, \R U — Rt\ < 1 < R u + R t , the number of Rt events is around 1%. 

Because the data are affected by errors we have to find all the possible forms 
for R u and R t , and for all the six unitarity triangles, and we can write a similar 
form for \\- Thus there is a similarity between the two approaches and in this 
form they have to lead to similar results. 

Let us see how the problem is solved in the usual manner [Hj, [Zj, [Hj, Hi, 

. Since the unitarity is assumed and built in the form Q the relation (|15|) is 
an identity; this means that we have to prescribe how the experimental data are 
related to the U entries. The first attempt was done by Wolfenstein ^S] , its nice 
feature being the possibility to estimate the order of magnitude for all U entries 
at a time when the experimental data were scarce. Being an approximation de- 
pending on a small parameter, the unitarity is far from being exactly satisfied. 
Taking into account that there is no simple relation between the parameters de- 
fined by Wolfenstein, A, A, p, ?y, and directly measurable quantities, Branco and 
Lavoura |19j proposed that instead of 5 it is better to choose another parameter, 
Q = \Utd/{U us U c b)\ 2 , this one being more relevant to B— mesons physics and a 
measurable quantity. Although the given parametrization is a Wolfcnstein-type 
one with four independent moduli almost nobody paid attention to it. 

Nowadays use of standard unitarity triangle approach is largely based on a 
paper by Buras et al |B] . They tried to restore the unitarity of the CKM matrix 
which in Wolfenstein parametrization is not satisfied exactly. Thus they choose 
four independent parameters as given by, see their Eq.(2.2), 

Sl2 = \Uus\, Sl3 = \U u b\, S 2 3 = \Ucb\, 6 (17) 
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and that set is replaced by 

A, A 7 p, 77 (18) 

for providing a unitary CKM matrix in terms of Wolfcnstein parameters. To 
this end they go back to H17J1 and impose the relations 

812 = A, s 23 =A A 2 , s 13 e- lS = A\ 3 (p-i v ) (19) 

to all orders in A. It follows that 

p = 513 cos 5, -q = Sl3 sin (20) 
S12S23 S12S23 

They observe that l|19f) and (|20[1 represent simply the change of variable from 
(I17fl to (|18|) . Further below they say: Making this change of variables in the 
standard parametrization we find the CKM matrix as a function of (A, A, p, 77), 
which satisfy unitarity exactly. So the unitarity problem of a parametrization a 
la Wolfenstein is completely solved. However . from the relations (|2()(l we get 



*12 6 23 V us v cb P 



The above relations tell us that only the first one is expressed in terms of 
measurable quantities, the second one is not. Taking into account that 6 is the 
parameter which describes directly the CP-violation it would been better to 
use for its determination a formula given in terms of measurable quantities. As 
we suggested at the beginning of this section this can be done, and a detailed 
comparison of both the methods, as well as its implications will be treated 
elsewhere [217] . 



4 Results 

For our fits in this paper we used the same 17 independent expressions for cos<5 
as in £Q , and for improving the statistics we considered another 7 groups of four 
independent moduli leading to a total of 42 independent expressions for cos 5. 
We considered both the experimental data ^2] and ^H], and the results have 
many features in common, although there are a few remarkable differences. 

Taking as test function \ 2 — Xi + x\ given as in Eas. p3() - (|14|l we obtained 
two (approximate) unitary matrices, one from the data |12j . and the second by 
using the novel data As in ^ we relaxed the condition (|14|l by taking all 
combinations with five and respectively four Vij that provided (^) + (5) =21 
new matrices. The set of these matrices was considered as 22 independent 
"experiments" for each experimental data set on which the statistics was done. 
For cos S these lead to 22 x 42 = 924 values that gave 6 = 89.979° ± 0.041° for 
Hagiwara et al data [T^, and to S = 89.9954° ± 0.00023° for Edelman et al data 
[T3*| . It is remarkable that the results of the fits concerning S are consistent with 
the Monte Carlo simulations, see FigQ] The \\ values are within the interval 
2. x 10 -7 < Xi — 1-7 x 10 -3 and the final results concerning the moduli are 
shown in the Table 1. 
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Quantity 


Central value dz error 


Central value zb error 




Data from [T3 


Data from 1131 




0.974868 ± 0.000031 


0.974534±0. 000027 


y US 


0.222717 ± 0.000031 


0.224138±0. 000027 




(5.356056 ± 0.035)10~ 3 


(3.34838±0. 005)10- 3 


v cd 


0.222529 ± 0.00003 


0.2240219±0. 000028 


Vcs 


0.9740447 ± 0.00005 


0.973699±0.000025 


v cb 


(4.1445 ±0.0064)10- 2 


(4.151237±0.0013)10- 2 


V td 


(10.60047 ±0.033)10- 3 


(9.862 ±0.0042)10~ 3 


Vts 


(4.042288 ±0.006)10~ 2 


(4.04626 ±0.0013)10- 2 


Vtd 


0.999126 ± 0.00008 


0.999132±0. 000014 



Tabic 1: Fit results and errors using the standard input from PDG data |12j - 
|13| . The errors are only statistical and come from the fact that x\ is n °t 
zero; for a unitary matrix the errors are zero. However they suggest the order 
of magnitude for the precision with which the experimental data have to be 
measured for a unambiguous determination of 8. The results show that there 
are unitary matrices compatible with the data. 

Looking at the final results we see that Vij values are not far from those 
provided by other fits, and are around the central values of both experimental 
data sets. The only unusual feature concerns the mean values for V ub , both 
being more or less far from the experimental values. In fact the intervals of 
variation were 1.9 x 10~ 4 < V u b < 3.72 x 10~ 3 , for Edelman et al. data [T3"| . 
and respectively, 3.29 x 10~ 3 < V u b < 11-6 x 10~ 3 , for Hagiwara et al. data 
|12j . The essential difference between the PDG data comes from the central 
values of the parameters showing that these ones cause two different behaviors 
for functions depending on Vij. But the most striking feature is the explicit 
strong dependence of the angle (3 on \U ub \. For seeing it we used the double 
stochasticity property and with the mean \/V matrices H and L given in Table 
1, obtained from Hagiwara et al. data, and, respectively, Edelman et al. data, 
we form the convex combination 

A = s/xH 2 + (1 - x)L 2 (22) 



matrix that depends on a continuous parameter. This dependence can be ex- 
pressed in terms of any Vij, at our choice. If we choose V u b as free parameter 
we find for sin 2/3 the behavior shown in Fig. 3. Taking the HFAG value ^1] for 
granted, sin 2/3 = 0.726 ± 0.037, we find 0.00380456 < V ub < 0.00440217 and 
the central value V ub = 0.00409197. Substituting them in Ea. f2l)l wc find three 
matrices, all of them being compatible to unitarity, the central one being 

/ 0.97464 0.223741 0.00409197 \ 
W= 0.22355 0.973809 0.0414911 (23) 
\ 0.0101019 0.04045 0.99913 / 

These results imply very strong correlations between all Vij and give small 
compatibility intervals for all Vij, e.g. we get 0.223542 < \U US \ < 0.223935, 
0.0414821 < \U cb \ < 0.0414998, 0.040447 < \U ts \ < 0.0404552, etc. Thus our 
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0.9 




0.0025 0.003 0.0035 0.004 0.0045 0.1 



Figure 3: sin 2/3 as function of V u b\ the horizontal band corresponds to sin 2/3 = 
0.726 ± 0.037 data, as given in [Hj. HFAG results imply 3.84456 x 10~ 3 < 
\U ub \ < 4.40217 x 10~ 3 . 

approach shows the constraining power of the unitarity method and allows us to 
obtain matrices compatible to unitarity and the experimental data in a direct 
way. 

5 Conclusions 

What we have done in the paper is a "numerical experiment" using the PDG 
data and the unitarity method for constraining the data, experiment that has 
shown a unexpected behavior of physical quantities, behavior that cannot be 
seen by using the previous methods for doing fits. As a conclusion we may say 
that our approach to the reconstruction of CKM unitary matrices being exact 
and using four independent moduli can be applied to both quark and lepton 
matrices and could outperform by far all the other methods used to reconstruct 
a CKM unitary matrix from experimental data. Preliminary results on the 
lepton matrix \/V |21j are very interesting and will be treated elsewhere. 

In the same time our approach seems to be the ideal tool for a coherent 
treating of all the experimental data available on moduli and angles of the 
unitarity triangles allowing to obtain results as those shown in Fig. 3. 

On the other hand the test x\ ~ 0; associated to physical values for cos 5, 
that expresses the full content of unitarity, has to be satisfied by all the present 
and future fits, and it have to be taken seriously into account by all the people 
working on unitarity triangle problems. 

All of us have to encourage a wider appreciation of different parameteriza- 
tions of unitary CKM matrices and the connections between them. Of course a 
selection of a specific set of angles and/or phases has no intrinsic theoretical sig- 
nificance because all the choices are mathematically equivalent; however a clever 
choice may shed more light on important qualitative or experimental issues. In 
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this sense we consider that we have to pay a special attention to the Aleksan et 
al. approach whose main point was: any phase-convention-independent 

product of CKM elements is a linear combination of such four independent 
phases, with integer coefficients, point which passed almost unnoticed. That 
means that all the fits have to provide values for these four fundamental fre- 
quencies because the used form of the CKM unitary matrix may depend on the 
physical process we want to measure, by choosing it such that the dependence 
of the decay amplitude on the CKM parameters should be as simple as possible. 
Having these fundamental frequencies it will be a simple way to translate the 
results obtained in a specific parameterization to results obtained in another pa- 
rameterization used by somebody else. For an other phase convention different 
from the usual one see for example Botella et al. |22| . 

We think that an important task for both experimenters and theorists is to 
improve all our experimental and theoretical tools we are working with, having 
in view the second-generation B-decay experiments at the LHC machine, and 
we consider our approach as a small step in this direction. 
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